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A new release of the Monte Carlo program Herwig-|--|- (version 2.5) is now available. This 
version comes with a number of improvements including: new next-to-leading order matrix 
elements, including weak boson pair production; a colour reconnection model; diffractive 
processes; additional models of physics beyond the Standard Model; new leading-order 
matrix elements for hadron-hadron and lepton-lepton collisions as well as photon-initiated 
processes. 
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1 Introduction 

The last major public version (2.3) of Herwig++, is described in great detail in [1-1]. As we 
did not produce a release note for version 2.4 we describe all changes since version 2.3 in this 
release note. The manual will be updated to reflect these changes and this release note is only 
intended to highlight these new features and the other minor changes made since version 2.3. 

Please refer to [1] and the present paper if using version 2.5 of the program. 

The main new features of this version are: new next-to-leading order (NLO) matrix ele- 
ments, including weak boson pair production; a colour reconnection model; diffractive processes; 
additional models of physics beyond the Standard Model; new leading-order matrix elements 
for hadron-hadron and lepton-lepton collisions as well as photon-initiated processes. 

In addition, the MCQNLO [o] program can now generate partonic configurations which can 
be showered and hadronized using Herwig+-|- to produce events with next-to-leading-order 
accuracy in the MCONLO approach to matching NLO matrix elements and the parton shower. 

1.1 Availability 

The new program version, together with other useful files and information, can be obtained 
from the following web site: 

http : //projects .hepf orge . org/herwig/ 

In order to improve our response to user queries, all problems and requests for user support 
should be reported via the bug tracker on our wiki. Requests for an account to submit tickets 
and modify the wiki should be sent to herwig@projects.hepforge.org. 

Herwig-|--|- is released under the GNU General Public License (GPL) version 2 and the 
MCnet guidelines for the distribution and usage of event generator software in an academic 
setting, which are distributed together with the source, and can also be obtained from 

http : //www.montecarlonet . org/ index .php?p=Publi cat ions/Guidelines 
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2 POWHEG 



The previous version of Herwig++ included for the first time a number of processes simulated 
according to the POWHEG NLO parton shower matching scheme [6, 7]. The current release 
builds on this internal library of NLO-accurate POWHEG simulations, adding doubly resonant 
W~^W~, W^Z^ and Z^Z^ pair production processes, fermionic Higgs decays and e~^e~ — )■ qq. 

This release also sees a major restructuring of the code for the parton shower, in particular 
the code handling the real emission in the POWHEG scheme, in order to make the imple- 
mentation of further processes in the POWHEG scheme and other extensions in the future 
easier. 

Previously, the hard and soft matrix element corrections were implemented in a set of ded- 
icated MECorrection classes inheriting from the MECorrectionBase class. Similarly, the real 
emission corrections in the POWHEG scheme were implemented in a set of HardGenerator 
classes inheriting from the HardestEmissionGenerator base class. While this had some advan- 
tages, in particular allowing the corrections to be applied regardless of how the hard processes 
was generated (for example the corrections could still be used with events read from a LHE file), 
it led to a significant replication of code between the HardGenerator and matrix element classes 
implementing the B function in the POWHEG scheme. We have therefore restructured the 
code so that both the matrix element corrections and real corrections in the POWHEG scheme 
are implemented in the matrix element or Decayer class implementing the hard scattering or 
decay process. In addition, the functionality of the separate PowhegEvolver has been merged 
into the base Evolver class. 

2.1 Vector Boson Pair Production 

The simulation of weak boson pair production employs the NLO computations of Refs. [8-10] 
and its construction has been documented in detail in Ref. [11], together with substantial 
comparisons to the MCFM and MC@NLO programs [12, 13]. These studies show excellent 
agreement with MCFM and generally good agreement with MCONLO. Two examples of these 
comparisons can be seen in Fig. 1. Note that spin correlations in the vector boson decays 
are included at tree level for all production modes and decay channels. The simulation of 
leptonically decaying vector bosons is enhanced, including leading-log photon emission effects 
resummed according to the YFS formalism [14, 15]. 

The NLO matrix element class, MEPP2VVPowheg, inherits from that of the leading-order 
matrix element, MEPP2VV, thus the selection of particular final states is performed using the 
Process interface. Any one of the five production modes W^W~ , W^Z^, Z^Z^, W~^Z^, and 
W~Z^, can be activated by setting Process equal to 1-5 respectively in the input files. Lastly 
we note that, as with all other Herwig++ POWHEG simulations, the default renormalization 
and factorization scale chosen for generating the underlying Born kinematics, i.e. the initial 
qq — >■ VV configuration, is given by the invariant mass of the colourless final-state system. In 
generating the hardest emission from this configuration, the renormalization and factorization 
scales used to evaluate the real cross section are set to the transverse momentum of the emission, 
as mandated by the POWHEG formalism. 

2.2 e+e- qq 

The next-to-leading order matrix element for e~^e" — > qq is implemented, including the masses 
of the quarks, using the POWHEG scheme in the MEee2gZ2qqPowheg class. The virtual correc- 
tions were taken from Ref. [16]. The real contribution to both the hardest emission and the B 
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(a) Cosine of the polar angle of a positron produced 
from a decaying Z*' boson, in its rest frame, in Z'^Z'^ 
pair production at the Tevatron (-^s = 1.96 TeV). 
The Herwig-|-+ POWHEG prediction is shown in red 
while the NLO and LO predictions of MCFM appear 
as black stars and blue dots respectively. 
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(b) Predictions for the boson spectrum in 

W'^W~ pair production at the LHC, assuming a 
hadronic centre-of-mass energy of y/s = 14 TeV. 
Leading-order Herwig++ results are seen in blue while 
those of MC@NLO and the Herwig++ POWHEG simu- 
lation are visible as black and dashed red lines. 



Figure 1: Comparison of Herwig-|--|-, MCFM and MCONLO for di-boson production. 

function was calculated using the internal helicity amplitude code. We use the massive dipole 
scheme of Ref. [17] to subtract the singularities from the real contribution for the calculation 
of B and to separate the singular regions for the generation of the hardest emission. 

The results for the cross section for both bottom and top production are in excellent agree- 
ment with the analytic results given in the appendix of Ref. [17]. The variation of the cross 
section for the production of bb pairs in e'^e~ collisions with the centre-of-mass energy is shown^ 
in Fig. 2(a). The results of this new simulation are in comparable, albeit marginally better, 
agreement with the data from the LEP experiments, for example the distribution of the thrust 
shown in Fig. 2(b), than the default Herwig-|--|- approach including a matrix element correction. 

2.3 Higgs Decay 

There are large QCD corrections to the partial widths for — )■ qq. As with the simulation of 
e~^e~ — )■ qq we use the dipole subtraction scheme of Ref. [17] both to subtract the singularities 
from the real contribution for the calculation of B and to separate the singular regions for 
the generation of the hardest emission. The internal helicity amplitude code was used to 
calculate the real emission and the results of Ref. [19] for the virtual contribution. In order 
to correspond as closely as possible with the masses used in the parton shower, we use the 
running mass for the coupling of the Higgs boson to the quarks, in order to resum the large 
corrections, and the pole mass in the kinematics. The resulting simulation is available in the 
SMHiggsFermionsPOWHEGDecayer class. 

2.4 Summary 

For convenience, we give the full complement of NLO-accurate POWHEG simulation classes 
included in Herwig-|--|-: 

^The analytic result is not shown as it is indistinguishable from the Herwig++ result. 
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Figure 2: Cross section for e'^e — )■ bb and thrust distribution at LEP. 



• the MEqq2gZ2ffPowheg and MEqq2W2ffPowheg classes for the production and decay of 
the 7*/Z° and bosons, respectively, in the Drell-Yan process; 

• the MEPP2HiggsPowheg class for the production of the Higgs boson via the gluon-gluon 
fusion process; 

• the MEPP2ZHPowheg and MEPP2WHPowheg classes for the production of the Higgs 
boson in association with the and bosons, respectively; 

• the MEPP2VVPowheg class for vector boson pair production processes; 

• the MEee2gZ2qqPowheg class for e^e~ — )■ qq; 

• the SMHiggsFermionsPOWHEGDecayer class for Higgs boson decay into quark-antiquark 
pairs. 

Examples illustrating the use of all of these POWHEG processes can be found in the 
LHC-Powheg. in and TVT-Powheg . in example files provided in the release. 



3 Colour Reconnection 

As of this release of Herwig-f-|-, a model for colour reconnection is included. It is implemented 
in the ColourReconnector class. The model can be regarded as an extension of the cluster 
model [20], which is used for hadronization in Herwig-|--|-. We would like to stress, however, 
that this colour reconnection model differs from the one used in (FORTRAN-) HERWIG [21], 
which is based on the spacetime structure of the event as described in Ref. [22]. 



3.1 Review of Cluster Hadronization 

Hadronization in Herwig-|--|- is based on the pre-confinement property of perturbative QCD [23]. 
According to that, a parton shower evolving to the cut-off scale Qq ends up in a state of 
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(a) Differential 2-jet rate, data from Ref. [24]. 
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Figure 3: Example distributions comparing an old tune of Herwig+-|- without colour reconnec- 
tions to a new tune including colour reconnections. 



colourless parton combinations with finite mass of 0{Qq). In the cluster hadronization model, 
these parton combinations - the clusters - are interpreted as highly excited pre-hadronic states. 
They act as a starting point for the generation of hadrons via cluster decays, which is possibly 
performed in multiple steps. This hadronization model is described in more detail in Ref. [1]. 

3.2 Colour Reconnection Algorithm 

Colour reconnection in the cluster model occurs at the stage where clusters are formed from 
the parton-shower products. Starting with the clusters that are produced generically by virtue 
of pre-confinement, the cluster creation procedure is slightly modified. This is done by allowing 
pairs of clusters to be 'reconnected'. That means the coloured constituent of cluster A and the 
anti-coloured constituent of cluster B form a new cluster, as do the remaining two partons. 
The following steps are performed for each cluster: 

1. loop over all other existing clusters and choose the one where a reconnection of the two 
clusters would result in the smallest sum of cluster masses; 

2. if such a reconnection possibility is found, accept it with probability 
ReconnectionProbability. 

3.3 Effects on Observables 

A re-tuning of the parton-shower and hadronization-related parameters to LEP data was in- 
evitable since this model for non-perturbative colour reconnection affects hadronization. We 
find agreement with the LEP data at a similar level to that achieved by the old Herwig+-|- 
tune without colour reconnection. As an example, the 2-jet rate using the Durham jet algo- 
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Figure 4: Comparison of Herwig-|--|- 2.4.2 and Herwig-|--|- 2.5 to ATLAS minimum-bias distri- 
butions at y/s = 0.9 TeV with iVeh > 6, > 500 MeV and \r]\ < 2.5. The ATLAS data were 
read off from plots pubhshed in Ref. [25]. 



rithm is shown in Fig. 3(a). A further illustration is given in Fig. 3(b), which shows the scaled 
momentum of charged particles at LEP. 

Since colour reconnection explicitly affects the cluster mass spectrum, the average multipli- 
city of charged particles changes. This observable was measured to be (Nch) = 20.92 ± 0.24 at 
LEP for ^/s = 91.2 GeV [18]. With (iVch)Hw++ = 20.73 Herwig++ (with colour reconnection) 
coincides with this measurement. 

Major improvements, however, are achieved in the description of minimum bias events and 
the underlying event in hadron collisions. In Fig. 4 we show the pseudorapidity distribution of 
charged particles and their average transverse momentum as function of the charged-particle 
multiplicity for minimum bias events at y/s = 0.9 TeV at the LHC This direct comparison of 
the current Herwig-|--|- release to its predecessor shows the progress in modelling minimum-bias 
events, partially enabled by the presence of non-perturbative colour reconnection. 

A second ingredient for the latest developments in the minimum bias and underlying event 
modelling is a modification of the model for soft interactions in Herwig-|--|- [26]. The previ- 
ous model generated additional soft scatters, implemented as soft gluon collisions in di-quark 
scatterings, in a way that their colour structure is entirely disconnected from the rest of the 
event. We extend this choice by allowing a colour connection between the soft scatters and 
the beam remnants. The parameter colourDisrupt is the probability of a soft scatter to be 
disconnected. 

Herwig-|--|- is also capable of reproducing minimum-bias data at 7 TeV, as well as underlying- 
event-related observables for both 0.9 TeV and 7 TeV. Different tunes for each application have 
to be used, though, since a proper description of all energies with the same set of parameters 
is not possible with the present model. A more generally applicable model is intended for 
future Herwig-|--|- releases. The default parameters in the current release remain those used in 
the previous version, i. e. no colour reconnections are included. The best parameters including 
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Figure 5: Diagram of a two-photon induced process. 



colour reconnections for different processes and energies can be found at 
http : //projects .hepf orge . org/herwig/trac/wiki/MB_UE_tunes 

4 Diffractive and Photon-Initiated Processes 

Some interesting forward processes in hadron-hadron collisions have been added in this version 
of Herwig++, namely photon-initiated processes and inclusive hard diffractive processes. 

4.1 Photon-Initiated Processes 

In photon-initiated processes photons are emitted by the incoming protons and their interaction 
yields a system X which is separated by large rapidity gaps in the forward region from the 
outgoing scattered protons. The protons remain intact in the interaction and they are only 
scattered by a small angle. A sketch of the interaction pp^p + 'y + j + p^p^X^p, where 
© stands for a rapidity gap, can be seen in Fig. 5. The emission of photons by protons is well 
described by Quantum Electrodynamics in the framework of Equivalent Photon Approximation 
[27]. The Photons are almost real (the photon virtuality ~ 0) and therefore the cross 
section can be factorized into the photon fluxes and the sub- matrix element. The photon flux 
is implemented in the BudnevPDF class. The photons are generated according to 



dN 



IT E^ Q 



(1) 



where E^ is the photon energy, E is the energy of proton and Q^in minimum virtuality 

of the photon allowed by the kinematics. The electric and magnetic form factors, Fe and Fm, 
are given by 

Fe = {4mlGl + Q^FM)/{4:ml + Q^), (2a) 
G| = FM/fil = il + QVQlr\ (2b) 

with the magnetic moment of the proton fip = 7.78, the fitted scale = 0.71 GeV^ and the 
proton mass rUp. An example input file example for LHC settings, LHC-GammaGamma . in, is 
provided in the release. 
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4.2 DifFractive Processes 

In the double pomeron exchange process both protons are left intact by the scattering process 
and two rapidity gaps occur in forward regions. In single diffractive processes one proton dissoci- 
ates, while the other is left intact, and a rapidity gap appears on one side of the event. Examples 
of both processes are shown in Fig. 6. The Ingelman-Schlein model [2<S] has been implemented 
to describe those processes. In this model the cross section is factorizes into a diffractive dis- 
tribution function and the cross section of the sub-process, a = /^(x, Q^, xp, t) ® (Tsuhix , Q'^) . 
The diffractive distribution function can be further decomposed into the pomeron/reggeon flux 
and distribution functions, respectively, 

/z)(x, xp, t) = fv/p(t, xp)fi/¥{Q^, x) + rj^fR/pit, XR)fi/R{Q^, x), (3) 

where /p(t, xp) and /R(t, xm) are pomeron and reggeon fluxes. We use the Regge theory inspired 
form 

/p(R)(t,a;p(M)) = Ap(M)e*«7x'("'««-"n«)*)-\ (4) 

where q;p(m)o and a'^-g^-^ are the pomeron/reggeon intercept and slope, respectively. Both the 
pomeron and reggeon fluxes are implemented in the PomeronFlux class because they differ only 
in choice of ao and a'. The pomeron distribution function fi/r>{Q'^,x) is implemented in the 
PomeronPDF class. Herwig++ includes the 2006 Fit A, 2006 Fit B and 2007 flts of the pomeron 
distribution function measured at HERA [29,30]. The PomeronPDF class allows the user to 
freeze or extrapolate the pomeron distribution functions outside the limits of the experimental 
fit. In the case of the reggeon distribution function, fi/M.{Q'^,x), an external (pion) parton 
distribution function from LHAPDF [31] can be used and Herwig++ provides only the interface 
class Reggeon PDF which mimics a reggeon beam particle and calls an external user-defined 
PDF. Interfaces to change parameters such as ap(K)o, ap(]R) and tjr are provided. By default 
they are set to the values determined in the used PDF fit. 

We also provide two options for the internal valence structure of the pomeron. The pomeron 
can be composed of either valence qq or valence gluons, the latter of which is the default setting. 
The reggeon is treated as a qq object since using the pion PDF is a usual and consistent approach 
to describe the reggeon. 

The pomeron or reggeon contributions can be simulated either separately or as a mixture. 
An example of the settings for the LHC can be found in the LHC-Dif f ractive . in input file 
supplied with the release. It should be noted that the model does not include the gap survival 
probability factors which need to be taken into account in hadron-hadron collisions. 

5 BSM Physics 

In addition to including more models of physics beyond the Standard Model (BSM), this release 
includes a number of other improvements to the simulation of BSM physics in Herwig++. 

5.1 Process Generation 

The inheritance structure of the code for simulating general or resonant 2—7-2 scat- 
tering processes has been changed. A new HardProcessConstructor base class has been 
added. The existing code for general 2 — )■ 2 scattering processes has been renamed 
TwoToTwoProcessConstructor and now inherits from the new HardProcessConstructor base class, 
as does the ResonantProcessConstructor class for resonant 2 — 2 scattering processes. As 
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Figure 6: Examples of single diffractive (left) and double pomeron exchange (right) processes. 

part of this change, the specific interfaces to the objects needed to simulate 2 — )■ 2 scat- 
tering processes have been removed from the ModelGenerator class and replaced by a list of 
HardProcessConstructors in order to make adding further types of hard processes easier. 

These changes have allowed the generation of hard scattering processes in BSM models to 
be extended to include some 2 — 3 processes involving neutral Higgs bosons. The production 
of Higgs bosons in association with an electroweak vector boson, and Z^, is generated using 
the HiggsVectorBosonProcessConstructor class. The production of a Higgs boson in association 
with a heavy quark-antiquark pair is calculated using the QQHiggsProcessConstructor class. The 
production of the Higgs boson via the vector boson fusion (VBF) process is simulated using 
the HiggsVBFProcessConstructor class. 

In addition, a number of improvements have been made to the hard process selection for 
2 — )■ 2 scattering processes in BSM models. The Processes interface can be used to se- 
lect between: [SingleParticlelnclusive] where at least one of the particles in the list of 
Outgoing particles must be produced, which is the previous behaviour and the new default; 
the [TwoParticlelnclusive] option where both the particles produced must be in the list of 
Outgoing particles; and the [Exclusive] option where only two particles are allowed in the 
list of Outgoing particles and both of these must be produced. In addition, the Excluded 
interface can be used to forbid specific particles as intermediate particles in the scattering pro- 
cess. Similarly the ExcludedVertices interface can be used to forbid specific vertices in the 
hard scattering processes 

The scale choice in the hard 2 — > 2 scattering processes in BSM models can now be changed 
via the new ScaleChoice interface to choose between s (default for colour-neutral intermedi- 
ates) and the transverse mass (default for all other processes). 

5.2 ADD Model 

In the ADD Model [32-38] gravity propagates in extra spatial dimensions, which have a fiat 
metric. The large size of these extra dimensions leads to a tower of Kaluza-Klein excitations of 
the graviton. They can either contribute as virtual particles to Standard Model (SM) processes, 
or be produced leading to missing energy signatures. This model is implemented using the 
conventions of Ref. [39]. 

An example input file for the LHC, LHC-ADD.in, is provided, and an example of the total 
cross section for photon or jet production with missing transverse energy is shown in Fig. 7. 
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(a) Photon + missing transverse energy cross section. (b) Jet + missing transverse energy cross section. 

Figure 7: The total photon and jet + missing transverse energy cross section at the LHC The 
photons/jets are required to have transverse energy greater than E'™™. The pseudorapidity of 
the photons/jets are required to satisfy, |?7^| < 2.5 and |?7jct| < 3.0, respectively. The cross 
section is integrated over either all values of s or s < M|,. This figure is based on Figs. 3 and 5 
ofRef. [39]. 

5.3 Leptoquarks 

Fermion masses may arise from the mixing of elementary fermions with composite, fermionic 
resonances of a strong sector [40] responsible for the breaking of the SU (2) l x U{1)y electroweak 
symmetry of the SM. It follows that this strongly-coupled sector must also be charged under 
colour SU{3) and must contain, at the very least, colour-triplet fermionic resonances that 
can mix with the elementary colour triplets to make the observed quarks. It is reasonable 
to expect that such a strongly-coupled sector will contain other coloured resonances. These 
may be bosonic and, depending on their charges, may couple to a quark and a lepton. These 
leptoquark resonances may be light if they arise as pseudo-Nambu Goldstone bosons and make 
an ideal target for LHC searches. They will decay exclusively to third-generation fermions due 
to suppression of the couplings to light fermions. 

The present implementation includes non-derivatively coupled leptoquarks, that couple to 
Standard Model fermions as in Eq. (2.4) of Ref. [11] and single- derivatively coupled leptoquarks 
such as those in Eq. (2.5) of Ref. [41]. In the case of the derivatively coupled leptoquarks, the 
simplification that the primed lepton and primed quark couplings are equal has been made. 

5.4 NMSSM 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) extends the Minimal Super- 
symmetric Standard Model (MSSM) with the addition of a singlet Higgs superfield, S. This 
leads to a larger particle content than the MSSM, i.e. 3 scalar Higgs bosons, 2 pseudoscalar 
Higgs bosons and 5 neutralinos. This model has been shown to overcome or reduce many 
problems associated with the MSSM. For recent reviews of the NMSSM see Refs. [42,43]. 

The phenomenologically relevant trilinear interactions of the Higgs bosons with fermions, 
sfermions, gauginos, gauge bosons, gluons, photons and the other Higgs bosons are included. 
All the other important interactions are inherited from the implementation of the Minimal 
Supersymmetric Standard Model. More details are available in Ref. [11]. 
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5.5 Transplanckian Scattering Model 

The interaction described by the Transplanckian matrix element is 2 — )■ 2 high-energy gravita- 
tional scattering of partons using the eikonal approximation [15]. The approximation is valid in 
the high-energy, low-angle scattering regime, where the centre-of-mass scattering angle of the 
incoming parton ^ — > or, in terms of Mandelstam variables, —t/s — t- 0. The implementation 
allows variation of the Planck scale, M^i, as well as the number of extra dimensions, up to a 
maximum of 6. 

5.6 Minor Changes to BSM Models 

The SusyBase class has been modified so that it can read the SLHA parameter file from the 
header of a Les Houches event (LHE) file. The name of the LHE file must be given via the 
setup command to the model class instead of the SLHA file. 

The vertices for the interactions of gravitinos in SUSY models have been added to allow 
the decay of the next-to-lightest supersymmetric particles (NLSP) in gauge-mediated SUSY 
breaking models. In addition, the approach used to calculate Rarita-Schwinger spinors for spin- 
I particles has been changed to improve the numerical stability for very light spin-| particles, 
such as the gravitino. 

Similarly, the vertex for the flavour-changing interaction of the stop quark, neutralino and 
charm quark has been added to allow the simulation of the decay of the lightest stop at pa- 
rameter points where its mass is lighter than that of the lightest chargino. 

In general, Herwig-|--|- does not include some three- and four-point vertices which do not 
have any phenomenologically relevant collider signals, for example three- and four-point self 
couplings of the Higgs boson and the four-point coupling of two vector and two Higgs bosons. 
Although these interactions are not relevant for the processes we include in the Standard Model, 
they can be in BSM models, due to the additional particle content. We have therefore extended 
the implementation of the Standard Model to include the triple Higgs boson self coupling and 
the coupling of two Higgs bosons to two vector bosons to make implementing these vertices in 
BSM models simpler. These vertices are also implemented in both the MSSM and NMSSM. 

The handling of the colour flows in 2 — 2 hard scatterings has been changed to improve 
the generality of the code and make future extensions easier. There should be no changes 
to the results for the currently implemented models. In addition, a number of improvements 
were made for the colour flows involving colour singlet exchange in the s-channel to allow the 
simulation of resonant Higgs boson production. 

The Helicity classes in ThePEG have been significantly cleaned up, leading to major sim- 
plification of the vertex implementations for all Standard Model and BSM physics models in 
Herwig++. The separate SpinBase and Spinlnfo classes have been merged into the new Spinlnfo 
class significantly simplifying the code structure and reducing the number of pointer casts. 
These changes, which should be transparent to most users, have been made to the helicity am- 
plitude code in order to improve the structure of the code and make future extensions easier. 

6 New Matrix Elements 

A number of new matrix elements are included in this release: 

• the MEPP2HiggsVBF for the production of the Higgs boson via electroweak vector boson 
fusion in hadron-hadron collisions^; 

^This matrix element is based on the dominant ^-channel contribution to this process. 



11 



• the MEPP2VV class for the production of pairs of electroweak gauge bosons, W^W , 
W^Z^ and Z^Z^ in hadron-hadron coUisions; 

• the MEPP2VGamma class for the production of an electroweak gauge boson, or Z^, 
in association with a photon; 

• the MEPP2QQHiggs class for the production of the Higgs boson in association with a tt 
or bb pair; 

• the MEee2VV class for the production of pairs of electroweak gauge bosons, W^W~ and 
Z'^Z^ in lepton-lepton collisions; 

• the MEGammaGamma2ff class for the production of fermion-antifermion pairs in photon- 
photon initiated processes; 

• the MEGammaGamma2WW class for the production of W^W~ pairs in photon-photon 
initiated processes; 

• the MEGammaP2Jets class for the production of jets in photon-hadron collisions via the 
partonic processes 57 — > qg^ (77 — > qg and g^ qq. 



7 Other Changes 

A number of changes have been made to the implementation of the SplittingFunctions in the 
parton shower module. Previously a splitting function class inheriting SplittingFunction had 
to be implemented for each combination of the colours and spins of the interacting particles. 
This has been changed so that now all the possible colour states of the interacting particles are 
implemented in the base class while the inheriting classes implement the spin structure only. 
This makes adding new interactions in BSM models simpler. In addition, the option of deleting 
specific types of branching in the shower has been added. 

A number of other more minor changes have been made. The following changes have been 
made to improve the physics simulation: 

• Soft QED radiation in Z^ decays is now fully 1-loop by default. As part of this change, 
the old WZDecayer class, which handled both and Z° decays, has been split into 
separate WDecayer and ZDecayer classes handling and Z^ decays, respectively. 

• Problems with numerical instabilities in the boosts applied to tau lepton decays have 
been fixed by postponing the tau decays until after the parton shower. A new interface 
setting Excluded in the DecayHandler class is available to prevent decays in the shower 
step. By default this is only enabled for tau leptons. 

• The default PDF set is now the LO* set of Ref. [46]. In addition, the PDFs used in the 
shower can now be set separately to those in the hard process using the ShowerHandler 
interface of the ShowerHandler class. 

• The shower, hadronization and underlying event parameters were retuned against LEP 
and Tevatron data respectively. 

• The mixing of B^ — B^ and B^ — B^ mesons has been added. The oscillation of the mesons 
is simulated including the CP-violating terms. However, there is no special treatment of 
those modes where both the meson and its antiparticle decay into the same final state. 
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• The missing colour structures required for a number of decays in BSM models have been 
added. 

• QED radiation is now enabled in all perturbative Standard Model and BSM decays and all 
perturbative decays by default, where there are no strongly interacting particles involved 
in the decay. 

• Spin correlations are now switched on by default for all perturbative decays. 

• New interfaces to the AcerDet [17] and PGS [48] fast detector simulations are now available 
in the Contrib directory. 

• FastJet [49] is now the only supported jet finder code. All analyses have been converted 
to use FastJet. 

• Improvements have been made to the momentum reshuffling in Deep Inelastic Scattering. 

• The LesHouchesReader class has been modified to allow processes that violate baryon 
number to be read from LHE files. 

• A number of changes have been made in the parton shower to allow the violation of 
baryon number in the hard process. The showering of baryon number violating decays 
was already supported. 

• Code for the simulation of the production of W'^ bosons using the POWHEG approach 
is now included in the Contrib directory. 

• fi~ , and z/e, and their antiparticles, are now available as beam particles. They are 
all supported in the DIS matrix elements, n'^^" collisions are supported in the general 
matrix element code for BSM models but not yet in the hard-coded matrix elements for 
lepton-lepton scattering. 

• The option of bottom and charm quarks is now supported for heavy quark production in 
the MEHeavyQuark class. 

• The polarization of tau leptons can now be forced in the TauDecayer class to assist in 
studies of tau polarization. 

A number of technical changes have been made: 

• The Tests directory has been added. It contains many additional input files to perform 
more detailed tests of the program than are performed by default during make check. 
These tests use both our own internal analyses and many of the analyses available in 
Rivet [oO]. 

• The KTgX output has been updated. After each run, a KT[t;X file is produced that contains 
the full list of citations. Please include the relevant ones in any publication. 

• A number of improvements for OS X systems have been made including fixes for the Snow 
Leopard release. 

• Makef ile-UserModules now includes the Herwig++ version number. In addition, the 
compiler flag -pedantic is no longer enabled so that user code using ROOT will compile. 
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• The obsolete KtJet and CLHEP interfaces have been removed. 

• The zero-momentum interacting particle used for bookkeeping in Min-Bias events is now 
labelled as a pomeron. 

• K^/K^ oscillations into ^ now occur at the production vertex of the kaon to give the 
correct decay length. 

• The default scale choice in POWHEG processes is now the mass of the colour- singlet 
system. 

• A new ZERO variable has been introduced, which can be used to set any dimensionful 
quantity to zero avoiding explicit constructs like . 0*GeV. 

• A number of improvements have been made to the implementation of three-body decays 
in BSM models. 

• The option of redirecting all output to stdout is now supported. The files previously 
ending in -UE.out and -BSMinfo.out are now appended to the log file. They now also 
obey the UseStdout flag. 

• The detection of FastJet in the configuration process has been improved. 

• The interfaces to Rivet and HepMC have been moved from Herwig++ to ThePEG. 

• The configuration process now looks for ThePEG in the location specified by — prefix. 

• Important configuration information is now listed at the end of the configuration process 
and in the file conf ig.thepeg. Please provide this file in any bug reports. 

• The Exception specifiers in the definition of the doinit() etc. member functions have been 
removed. This may require the removal of the exception specifier after the function name 
in some user code. 

• The SLHA EXTPAR block can now be used to set tan/3. 

• The warning threshold for branching ratios not summing to 1 has been relaxed. It is now 
a user interface parameter. 

• The cross section for inelastic scattering for minimum bias processes is now available in 
the -UE . out files. 

• A number of classes have been renamed so that tilde is correctly spelt. 

• Some deprecated interfaces in the MPIHandler class have been removed. 

• The unused doubly heavy baryons have been deleted from the input files. 

• A number of minor changes have been made to improve the stability of the phase-space 
integration in particle decays. 

• The handling of the debugging flags has been made more consistent. The user now needs 
to add, for example, set LHCGenerator :DebugLevel 1 in the input file or the -d 1 
option on the command line in order for any debugging printout, including printing of 
events to the log file, to be enabled. 
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• The handling of the PDG codes for particles has been changed in order to prevent prob- 
lems on 64-bit systems mapping unsigned int to long. 

• The version of libtool has been updated to 2.4. We also use the silent rules available in 
recent versions of automake, reducing the output during the compilation of Herwig+-(-. 

• The version of LoopTools used in Herwig++ has been updated to 2.6. 

• A number of .ice files have been removed and the corresponding code moved to either 
the corresponding .h or . cc files. Similarly a number of unnecessary inline directives 
and unused header files have been removed. 

• The implementation of the HerwigRun class has been simplified by using more features of 
the EventGenerator class of ThePEG. 

• A number of changes have been made to fix warning messages from the Intel compiler. 

• The inheritance structure of the SMHiggsWidthGenerator has been improved to bring it 
into line with other WidthGenerator classes. 

• The option of applying a sequential longitudinal and then transverse boost in the 
QTildeReconstructor in order to conserve energy and momentum, rather than one sin- 
gle boost, has been introduced primilarly for use by the MC@NLO program. 

• The interface to the LHAPDF package [31] is now a separate dynamically loadable module 
ThePEGLHAPDF . so rather than being part of libThePEG . so. 

• The irrelevant three-body decays are no longer included by default when using the 
Randall- Sundrum model. 

• The LeptonLeptonRemnant class has been renamed UnResolvedRemnant as it is now used 
in other applications, for example in photon radiation from protons in hadronic collisions. 

• The handling of weighted events in the internal AnalysisHandler classes has been improved. 

• Various changes have been made to allow Interpolator objects to be written to a persistent 
ThePEG stream. 

• The obsolete BSMCascadeAnalysis has been removed. 

• A value of the Fermi constant, Gp-, is now available in the StandardModelBase class. This 
allowed the removal of several local values of this parameter. 

• The handling of initial-state radiation when the incoming hadron is not along the 2;-axis 
has been improved. 

• The tolerance parameter used to check momentum conservation in the BasicConsistency 
analysis class can now be changed by the user. 

• A new switch has been added to the ModelGenerator class so that the branching ratios it 
calculates in BSM models can be outputted in the SLHA format. 



15 



The following bugs have been fixed: 

• The example input files for Powheg processes now set the NLO as correctly, and are run 
as part of make check. 

• A problem that led to the truncated shower not being applied in some cases has been 
fixed. 

• An accidental duplication in the calculation of event shapes was removed, they are now 
only calculated once per event. Several other minor issues in the event shape calculations 
have also been fixed. 

• An initialization problem in the internal MRST PDFs was fixed. 

• The scale in the Vertex classes can now be zero where physically possible. 

• The Herwig++ main program now correctly treats the -N fiag as optional. 

• The accuracy of boosts in the ^-direction has been improved to fix problems with ex- 
tremely high pt partons. 

• A bug in the implementation of the PDF weight in initial-state q qg splittings has 
been fixed. 

• A bug in the x^x^W^ and charged Higgs-sfermions vertices has been fixed. 

• The longitudinal boost of the centre-of-mass frame in hadronic collisions is correctly 
accounted for now in the generation of QED radiation. 

• Numerical problems have been fixed, which appeared in the rare case that the three- 
momenta of the decay products in two-body decays are zero in the rest frame of the 
decay particle. 

• The numerical stability in the RunningMass and QTildeReconstructor classes has been 
improved. 

• The stability of the boosts in the SOPTHY code for the simulation of QED radiation has 
been improved. 

• A problem with forced splittings in the Remnant was fixed. 

• A problem with emissions from antiquarks in the soft matrix element correction in 
e'^e~ —!> qq was fixed. This also required a new tune of the shower and hadronization 
parameters. 

• The matrix element correction for QCD radiation in decays, which was not being 
applied, is now correctly used in decays. 

• The presence of top quark decay modes in SLHA files is now handled correctly. 

• Additional protection against problems due to the shower reconstruction leading to par- 
tons with X > 1 has been added. 

• Changes have been made to allow arbitrary ordering of the outgoing particles in BSM 
processes. 
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• Two bugs involving tau decays have been fixed. The wrong masses were used in the 
KPiCurrent for the scalar form factors and a mistake in the selection of decay products 
lead to — 7- Tr^K~ being generated instead of — ?■ riK~ . 

• To avoid crashes, better protection has been introduced for the case where diquarks cannot 
be formed from the quarks in a baryon-number violating process. In addition, the parents 
of the baryon-number violating clusters have been changed to avoid problems with the 
conversion of the events to HepMC. 

• A bug in the QEDRadiationHandler class that resulted in no QED radiation being generated 
in W~ decays has been fixed. 

• A number of minor fixes to the SUSY models have been made. 

• A fix for the direction of the incoming particle in the calculation of two-body partial 
widths in BSM models has been made. 

• The LoopTools cache is now cleared more frequently to reduce the amount of memory 
used by the program. 

• Negative gluino masses are now correctly handled. 

• A problem with mixing matrices that are not square has been fixed. 

• A problem in the matrix element correction in e^e~ — )■ ti events has been fixed. 

• The MEee2gZ2ll has been fixed to only include the photon exchange diagram once rather 
than twice as previously. 

• A problem has been fixed that occurred if the same particle was included in the list of 
DecayParticles. 

• A number of minor problems in the vertices for the UED model have been fixed. 

• The missing identical-particle symmetry factor in MEPP2GammaGamma has been in- 
cluded. 

• A floating point problem in the matrix element correction for top decays has been fixed. 



8 Summary 

Herwig++ 2.5 is the seventh version of the Herwig-|--|- program with a complete simulation 
of hadron-hadron physics and contains a number of important improvements with respect to 
the previous version. The program has been extensively tested against a large number of 
observables from LEP, Tevatron and B factories. All the features needed for realistic studies 
for hadron-hadron collisions are now present and we look forward to feedback and input from 
users, especially from the Tevatron and LHC experiments. 

Our next major milestone is the release of version 3.0, which will be at least as complete 
as HERWIG in all aspects of LHC and linear-collider simulation. Following the release of 
Herwig++3.0, we expect that support for the FORTRAN program will cease. 
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